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Recently, D. B. M. Dickerscheid et al. have proposed 
a hamiltonian to describe bosons in an optical lattice 
across Feshbach resonance A similar hamiltonian 
has also been used by some of these authors later on to 
describe resonantly interacting Bose-Fermi mixtures Q. 
These hamiltonians are strictly proposals and are not 
based on systematic derivations. Here, we would like to 
point out that there are serious logical inconsistencies in 
these proposals as well as in the calculations performed 
on them. As a result, these models 0,0 can not describe 
resonantly-interacting quantum gases in optical lattices, 
and the calculations in these works are incorrect even if 
the errors in the construction of the model are ignored. 

The Hamiltonian in ref.[l| is H = hi+T+h c , where 
T is the tunneling between wells, h c is a sum of chemical 
potential terms, and hi is the single site Hamiltonian 
which is of the form (with the site index i suppressed 
and ignoring the background interactions irrelevant for 
our discussions), 

h = e a a)a + £ a b\b a + g' \[Z~l{b\aa + h.c). (1) 
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where a is a boson in the open channel restricted to the 
lowest band e a ; e a are the energies of the lowest two 
"dressed" states b a in a trap, cr =T, J,, which are linear 
combinations of the closed channel molecule d and a pair 
of bosons a 

bl = ^fZrd) T v 7 ! - Z~«cJa\ (2) 

where Z\ = 1 — Zj. if one restricts to a single band. The 
authors also stated explicitly (at the end of the second 
paragraph of Section II in ref. [l| ) that the dressed states 
are the two lowest states of a two-particle system in a 
single harmonic well, meaning that the Hamiltonian for 
a single well with only two particles is 

hsingle = 'Y] Z<jb\b a (3) 
a 

and nothing else. The problem of eq.Q) is that it can 
never reduce to the correct Hamiltonian h S i ng ie when it 


is reduced to two-particle in a single well, for g' and e a 
are non-zero. 

There has never been a real derivation of how the terms 
g' comes about in As it stands, it describes the con- 
version of the a pair of particles (aa) to one of the two 
dressed states. This is, however, not how a Feshbach res- 
onance works. Feshbach resonance is caused by the reso- 
nance between a pair of particles aa in the open channel 
with a tightly bound pair (or molecule) in closed channel 
d, and is described by the well known resonance model 

h res = e a a) a + vd^ d + g'(d^aa + h.c.), (4) 

where v is the detuning between particle pairs in the 
closed and open channel. The two-particle eigenstates of 
eq.ljl)! are precisely bps in eq.J2J) with Z a given by some 
function of v, e a and g' . Although the g' term in eq.(0} 
is the same g' term in eq. Q since d = ^fZ^b^ + ^fZb^ , 
eq. Q is still h s i ng i e after a and d are expressed in terms 
of b a . Writing h res as h amounts to double counting 
many terms, and is the reason why it cannot reduce to 
the correct two-body Hamiltonian. 

The same logical inconsistency in ref.(Q]) extends to 
its method of solution. The authors have ignored the 
relation (cq.(0)) when solving the lattice model. This 
amounts to treating the three states 6j, by, and a 2 on 
each site as independent two-particle states. This can 
not be because the dimensionality of the Hilbert space 
for the lowest doublet of two-particle in a well is precisely 
two. The failure to implement the relation eq.Q creates 
an unphysical Hilbert space on every site, thereby intro- 
ducing highly uncontrollable errors in the calculation. 

The fact that the authors have double counted many 
terms in the correct Hamiltonian, and have completely 
ignored the relations between dressed states and open 
channel states have led us to the conclusion stated in the 
first paragraph. The same errors also repeat in ref. 2}. 
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